Exact non-Markovian cavity dynamics strongly coupled to a reservoir 
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The exact non-Markovian dynamics of a microcavity strongly coupled to a general reservoir at 
arbitrary temperature is studied. With the exact master equation for the reduced density operator 
of the cavity system, we analytically solve the time evolution of the cavity state and the associated 
physical observables. We show that the non-Markovian dynamics is completely determined by the 
propagating (retarded) and correlation Green functions. Compare the non-Markovian behavior at 
finite temperature with those at zero-temperature limit or Born-Markov limit, we find that the non- 
Markovian memory effect can dramatically change the coherent and thermal dynamics of the cavity. 
We also numerically study the dissipation dynamics of the cavity through the mean mode amplitude 
decay and the average photon number decay in the microwave regime. It is shown that the strong 
coupling between the cavity and the reservoir results in a long-time dissipationless evolution to the 
cavity field amplitude, and its noise dynamics undergoes a critical transition from the weak to strong 
coupling due to the non-Markovian memory effect. 

PACS numbers: 42.50.-p, 03.65.Yz, 42.79.Gn 



I. INTRODUCTION 



The dissipation quantum dynamics of optical cavities 
has been well investigated and deeply understood un- 
der the Born-Markov (BM) approximation [1[. The BM 
approximation is valid when the coupling between the 
system and the environment is weak enough so that the 
perturbation is applied, and meantime the characteris- 
tic time of the environment is sufficiently shorter than 
that of the system so that the non-Markovian memory 
effect is negligible. However, in many situations in the 
recent development of optical microcavities, the strong 
coupling effect or the long-time memory effect has be- 
come an important factor in controlling cavity dynam- 
ics. Typical examples include optical fields propagating 
in cavity arrays or in an optical fiber |3-|J1 , trapped ions 
subjected to artificial colored noise [5|i8|, and microcav- 
ities interacting with a coupled resonator optical waveg- 
uide (CROW) or photonic crystals [9l-[la|. etc. Specif- 
ically, for the trapped ions coupled with an engineered 
reservoir, the change of the characteristic frequency of 
the reservoir can be accomplished simply by applying 
a random electric field through a band-pass filter defin- 
ing the frequency spectrum of the reservoir [5|. While, 
for a cavity interacting with CROW or photonic crys- 
tals, the coupling between them is controllable by chang- 
ing the geometrical parameters of the defect cavity and 
the distance between the cavity and the CROW p^ . 
Both of them provide non-Markovian dissipation and de- 
coherence channels [Q-ll, [l^ ■ These strong coupling or 
long-time memory effects result in a complicated non- 
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Markovian process in cavity systems that has become a 
crucial concern for the rapid development of quantum 
information and quantum computation in terms of pho- 
tons [ItI- The non-Markovian behavior of the trapped 
ions has been discussed in many works p-Q- ^^ this 
paper, we shall investigate the non-Markovian dynamics 
for the second case, the cavity strongly coupled with its 
environment (CROW or photonic crystals). 

Quantum dynamics of cavity systems is completely de- 
scribed by the master equation of the reduced density op- 
erator by taking the cavity as an open system. The mas- 
ter equation under the BM approximation can be found 
in many textbooks [l|, [iM[3- However, an exact master 
equation beyond the BM approximation is only derived in 
a few works. The first exact master equation, also called 
as Hu-Paz-Zhang master equation for quantum Brown- 
ian motion, was found almost two decades ago [20|, in 
terms of Wigner distribution function in_ phase space via 
Feynman- Vernon influence functional |21i] . Recently, our 
group developed exact master equations of the reduced 
density operator for both the fermion systems [22, [2j| 
and the bosonic systems ^2^ l25J, by using an extended 
Feynman- Vernon influence functional. These exact mas- 
ter equations can fully depict the quantum dissipative 
and decoherence dynamics in various open systems in the 
strong coupling regime. For the cavity system, the mas- 
ter equation obtained in [2J, |23 is for zero-temperature. 
Here we shall extend it to a finite temperature. In prac- 
tical, the temperature effect is also unavoidable and non- 
negligible. It has been pointed out ilj] that for cavity 
frequency lies in the microwave regime, thermal photons 
are presented even at liquid helium temperatures \2(S\. 
In this paper, we shall use the exact master equation 
which is valid at arbitrary temperature to investigate 
the exact non-Markovian dynamics of a general cavity 



strongly coupled with its reservoir, and to find general 
features of the coupling and temperature dependence in 
non-Markovian dynamics. 

The paper is organized as follows. In Sec. II, we in- 
troduce the exact master equation we developed recently 
for the reduced density operator of the cavity system cou- 
pled to a general reservoir, from which the second-order 
perturbation master equation and the BM maser equa- 
tion are reproduced at well-defined limits. In Sec. Ill, 
we then solve analytically the exact master equation via 
the coherent-state representation, included the exact so- 
lutions of mean field amplitude and the average photon 
number inside the cavity that experimentally measur- 
able. We also obtain explicitly the reduced density ma- 
trix for three different initial states: the vacuum state, 
the coherent state and the mixed state, to show the dif- 
ferent non-Markovian behaviors. The decoherence dy- 
namics of cavity field is explicitly analyzed. In Sec. IV, 
we numerically demonstrate the exact non-Markovian be- 
havior through the time-dependence of the mean mode 
amplitude and the average photon number in the cav- 
ity in both the weak and strong coupling regimes for 
three typical spectral densities, the Ohmic, sub-Ohmic 
and super-Ohmic cases, with the cavity's frequency be- 
ing focused in the microwave regime. We find that the 
strong coupling between the cavity and the reservoir re- 
sults in a long-time dissipationless evolution to the cavity 
field amplitude, and the noise dynamics undergoes a crit- 
ical transition from the weak to strong coupling due to 
the non-Markovian memory effect. Finally a conclusion 
is given in Sec. V. 



A. Exact master equation 

The exact master equation for the cavity field is given 
in terms of the reduced density operator which is de- 
fined from the density operator of the total system by 
tracing over entirely the environmental degrees of free- 
dom: p{t) = trRptot('^), where the total density operator 
is governed by the quantum Liouville equation ptot(i) = 
^-j-H(t-to) p^^^^^^yj■H(t-to) ^ As usual [13|, assuming 
that the cavity field is uncorrelated with the reservoir 
before the initial time Iq: ptot{to) = pita) (8> Pnito), 
and the reservoir is initially in the equilibrium state: 
Pfi(^o) = ■i-e"^^/'^'^^"^^ where ks is the Boltzmann con- 
stant and T the reservoir's initial temperature. Then 
tracing over all the environmental degrees of freedom 
can be easily carried out using the Feynman- Vernon in- 
fluence functional approach [21] in the framework of co- 
herent state path- integral representation [2^. The re- 
sulting master equation for the reduced density operator 
has a standard form similar to the master equation for 
electrons in nanostructure we developed recently [23, ^^ 
(with some sign difference due to the different statistical 
property between fermions and bosons): 

pit)^-iu;'oit)[a^a,p{t)] 
+ K (t) {2ap{t)a^ - a)ap{t) - p{t)a) a) 
+ K{t){a^p{t)a + ap{t)a^ - a) ap{t) - p{t)aa)], (2) 

where the time-dependent coefficient a;g(i) is the renor- 
malized frequency of the cavity, while K[t) and K(i) de- 
scribe the dissipation and noise to the cavity field due 
to the coupling with the reservoir. These coefficients are 
non-perturbatively determined by the following relations: 



II. EXACT MASTER EQUATION FOR A 
CAVITY IN A GENERAL RESERVOIR 



We consider a cavity with a single mode coupled to a 
general reservoir. The Hamiltonian of the total system is 
given by 



uj'^{t) = -lm[u{t)u-^ {t)], (3a) 

K(t)^-Re[u{t)u-^{t)], (3b) 

K(t) = v{t) - 2w(t)Re[u(t)u~^(i)], (3c) 

and u{t) and v{t) satisfies the integrodifferential equa- 
tions of motion: 



H = hwQa* a 



E 

k 



fujkblbk 



J2f^yk[bla + a^bkj, (1) 

k 



in which the first term is for the single cavity mode with 
a ' , a being the creation and annihilation operators of the 
cavity field, luq is its frequency; the second term is the 
Hamiltonian Hj^ for a general reservoir modeled as a col- 
lection of infinite harmonic oscillators, where bj, and b^ 
are the corresponding creation and annihilation operators 
of the kth oscillator with the frequency ujk ■ The coupling 
between the cavity and the environment is described by 
the third term and Vk is the coupling strength between 
them, which is tunable for the reservoir being CROW or 
photonic crystals J9|-KL6|. 



u{t) + iuJo'U'{T)+ / dr'gir - t')u{t') = 0, 



(4a) 



v{t) + iijJov{T)+ / dr' g{T — t')v{t') 



[ dT'g{T-r')u*{T), (4b) 

Jtn 



subjected to the initial condition u(to) = 1, and u(t) = 
u{t -I- io — t). Note that the integral kernels in the above 
equations involve non-perturbatively the time correla- 
tion functions of the reservoirs: g{T — r') and ^(t — r'). 
These two time-correlation functions characterize all the 
non-Markovian memory structures between the cavity 
and the reservoir. By defining the spectral density of 



the reservoir J{oj) — Svr^j, |Vfe|^(5(a; — Wfc), the time- 
correlation functions can be expressed as 



Jo ^'^ 



(5a) 



did 



9ir-r')= I ^JMn(^,T)e— (---'), (5b) 

Jo ^7^ 

where ri{uj,T) — nu,/k'gT_^ is the average number dis- 
tribution of the reservoir thermal excitation at the ini- 
tial time to- If the reservoir spectrum is continuous, 
Vk -^ ^(w), we have J(cl') — 2TTg{uj)\V{uj)\'^ where g{uj) 
is the density of state for the reservoir. 

The master equation ([2]) is exact, far beyond the BM 
approximation widely used in quantum optics. The back- 
reaction effect between the system and environment is 
fully taken into account by the time-dependent coeffi- 
cients, Wo(i), K(t) and K(t), in the master equation ([2]) 
through the integrodifferential equations (|H). In fact, we 
can directly solve from Eq. ^ 

u{t)u^^ {t) — —iujQ — I dTg{t ~ T)u{T)u^^{t), (6a) 

Jto 

v{t)= I rfri / (iT2 u(Ti)g(Ti -T2)u* (t2) . (6b) 
Jto Jto 

Then these time-dependent coefficients in the exact mas- 
ter equation can be simplified as 

uj'f^{t)=ujn+lm[w[t)], K{t)^Re[w{t)], (7a) 

K{t) = v{t) + 2v{t)Re[w{t)]. (7b) 

with w{t) = J^ dTg{t — t)u{t)u~^ (t) , which can be calcu- 
lated by solving u(r) non-perturbatively from Eq. dJa]). 
Thus, the non-Markovian memory structure is non- 
perturbatively built into the integral kernels in these 
equations. The expression of the integrodifferential equa- 
tion ((H) shows that u{t) is just the propagating func- 
tion of the cavity field (the retarded Green function in 
nonequilibrium Green function theory [29[), and v{t) is 
the corresponding correlation (Green) function which is 
also determined by u{t) [see the equation (j6b|) ]. There- 
fore, the exact master equation for cavity reduced density 
operator depicts the full nonequilibrium dynamics of the 
cavity system. 

In fact, the exact master equation presented here sim- 
ply covers the exact master equation at zero-temperature 
we derived very recently [24 . |25|. Taking the zero- 
temperature limit r = 0, then n{uj,T) = so that 
g{T — r') = 0. As a result, we have v{t) =0 and therefore 
the coefficient Ti{t) = 0. The master equation is simply 
reduced to 

+ K (i) {2a/9(i)at - a''ap{t) - p{t)a''a}, (8) 

with the temperature- independent coefficients WQ(i), K{t) 
obeying the same equation paj|3b[) through Eq. (Hal) . 



Eq. ([8]) is the exact master equation for cavity fields cou- 
pled to the vacuum fluctuation (i.e. the zero temperature 

limit) mn. 



B. Reduce to BM limit 

Interestingly, the exact master equation ([2]) has the 
same form as the BM master equation in the literature 
jT]. The difference is the coefficients in the master equa- 
tion. Here all the coefficients are time-dependent and 
determined by integrodifferential equations of motion (|4]) 
or ([6|), which non-perturbatively takes into account the 
back-reaction effects of the reservoir on the cavity field. 
While the coefficients in the BM master equation are all 
time-independent that ignores all the memory effects be- 
tween the cavity and the reservoir and are determined 
under the perturbation approximation up to the second 
order of the coupling V{ll!) and then taking the Markov 
limit. 

Explicitly, since the time correlation functions g{T — T') 
and g{T — r') are already proportional to |V^(w)p, taking 
approximately the time-dependent coefficients in Eq. ([7]) 
up to the second-order of the coupling V{lo) means that 
the propagating functions u{t) and u~^{t) in the right 
hand side of Eq. ([6|) should be approximated only up 
to the zero-order: uo{t) = ^-^'^oir-to) ^^^^ UQ^{t) = 
gjc^o(t-to). This leads to 



t poo 



u{t)u-^{t)^-iLJo- I dT I ^J(t^)e^''("-"«'(*-"), 
Jto Jo 27r 

(9a) 



/"* f°" duj 

v{t)^2 dr ■—J{uj)n{uj,T)cos[{uj-ujo){t-T)]. 

Jto Jo 27r 

(9b) 

The term 2v{t)Re[w{t)] in JTbl is proportional to \V{cv)\^ 
and should be ignored in the same approximation. Then 
the coefffcients of Eq. ([3]) or ([7]) are reduced to 

/"* 1"°° dui 
U!Q{t)~uJo- dr --J{uj)fim[{uj-uJo){t-T)], 

Jto Jo 27r 

(10a) 

/"* f°° du) 
K{t)~ dr — J(w)cos[(w-wo)(i-r)], (10b) 

Jto Jo 27r 

K(i) ~ 2 / dr — J{uj)n{uj,T)coa[{uj-ujo){t-T)]. 

Jto Jo 27r 

(10c) 

Substituting these coefficients into Eq. ([2]) results in the 
master equation of the cavity field in the perturbation 
approximation up to the second order of the coupling 
constant between the cavity and the reservoir. This per- 
turbative master equation is a good approximation only 
for the dissipation and noise dynamics of the cavity mode 
in the weak coupling regime. 



To reproduce the standard BM master equation with 
the time- independent coefficients, one needs further to 
take the Markov hmit where t is the typical time scale 
for the dynamics of the system and the t' integration is 
dominated by much shorter time characterizing the decay 
of reservoir correlations [i|, [23 • In other words, one can 
take t' integration to infinity in the equation pi 



lim 



UJ-UJq' 



(11) 



where V denotes the principle value of the integral. As 
a result, all the coefficients in the master equation, given 
by (|10p . become time- independent: 

(jj'q — ujQ + Soja, K — TTg{uja)\V{uja)\ = J{uj())/2, (12a) 
K = 2TTg{wo)\V{uo)fn{uQ,T) = 2Kn{uQ,T), (12b) 



°o J, .3(")|V(")|- 



Then 



and the frequency shift Sujq = V L duj 
the exact master equation ([2]) is reduced to 

p{t) = --(wo + SujQ)[a^ a, p{t)] 

+ K\2ap{t)a* — a*ap{t) — p{t)a'a\ 

+ 2Kn{u!o, T) \a* p{t)a + ap{t)a* — a* ap{t) — p(t)a(^\ . 

(13) 



Eq. p2|) with (fT3l) reproduces exactly the BM master 
equation in quantum optics for a single cavity mode in- 
teracting with a thermal field [l|. 



III. EXACT SOLUTION OF THE MASTER 
EQUATION 

In this section, we shall present the exact solution of 
master equation for some physical observables and also 
for the reduced density operator of the cavity. 



A. Exact solution for some physical observables 

The main physical observables for a cavity are the de- 
cays of the mean mode amplitude and the average photon 
number inside the cavity. The mean mode amplitude of 
the cavity field is defined by (a(i)) =tr[ap(i)]. From the 
exact master equation ([5]), it is easy to find that (a(i)) 
obeys the equation of motion 

{a{t)) = -KW + <mait)) = ^(« (t))- (14) 



which has the exact solution: 



{a{t))=u{t){a{to)). 



(15) 



That is, the time evolution and the decay behavior of 
the mean mode amplitude is totally determined by u{t). 



Eq. p^ clearly indicates that u(t) is the propagating 
function characterizing the time evolution of the cavity 
field. 

Another important physical observable is the average 
particle number inside the cavity, which is defined by 
n{t) =tT[d''ap{t)]. From the exact master equation, it is 
also easy to find that 

n{t) = -2K{t)n{t)+Ti{t). (16) 

On the other hand, Eq. pb[) can be rewritten as 



V{t) = ~2K{t)v{t)+K{t), 



(17) 



with —2K{t) = \u/u{t) + H.c.]. Combing these equations 
together, we obtain the exact solution of n{t) in terms of 
u{t) and v(t): 



n{t) ^u{t)n{to)u*{t)+v{t). 



(18) 



This relationship is similar to the fermion case we derived 
recently [30[. In fact, the above solution is a result of the 
correlated Green function in nonequilibrium Green func- 
tion theory [20|. Since v{t) is also determined by u{t) 
as we can see from Eq. (j6bl) . both the mean mode am- 
plitude and the average photon number inside the cavity 
are completely solved by the propagating function u{t). 
If we take the BM limit in which all the coefficients in 
the master equation are reduced to cjq = wo + Sluq , k = 
J{ujq)/2,k = 2i<n{ujQ,T) with Jwq — P Jq 



00 du Jji^) 
27r uj—uJo 



see 



Eq. ([T^ ]. then the resolution of u (t) simply becomes 



UBM (i) = e-(''"«+'')(*~*''^ 



(19) 



Correspondingly, the evolution of the mean mode ampli- 
tude in the BM limit is 



{a{t)) 



BM 



-{iu:'g+K){t~to) 



a (to)), 



(20) 



which shows an exponential decay in time. Substituting 
(|T9| into (gbl, we have the BM solution of v{t): 



z;BM(i)=^(^o,r)[l-e-2«(*-*o)] 



(21) 



Thus the average photon number in the BM limit is sim- 
ply given by, 

UBuit) = n(to)e-''=(*"*'') + n (c^o, T) [l - e-2K(*-t„)] 

(22) 
These results reproduce all the BM solutions in weak 
coupling regime [l]. However, the exact solutions of 
Eqs. dTSl) and (ITS)) allow us to explore the non-Markovian 
dynamics of the cavity systems not only in the weak 
coupling regime but also in the strong coupling regime. 
The explicit difference between the non-Markovian and 
Markov dynamics can be seen by comparing the solution 
of Eqs. (US]) and ^ with ^ and 



TABLE I: CoefBcients in the exact solution of the propagating function 





Ait) 


B{t) 


C(t) 


D{t) 


Exact result 


1 


l+v(t) 


v(t) 
l + v(t) 


1 l"(t)l^ 

l + v{t) 


l+v(t) 


T ^0 hmit 


1 


u{t) 





l-\u{t)f 


BM limit* 


1 


"bmW 

l + v-BM(t) 


i'Bm(*) 
l+"BM(t) 


1 I«bm{*)I^ 

l + "BM(t) 


l + "BM(t) 



* where UBuit) and VBuit) are given by Eqs. (|19p and ((21 



B. Exact Solution of the reduced density operator 

In fact, the reduced density operator can be also ex- 
plicitly obtained through the coherent state representa- 
tion. The reduced density matrix in the coherent state 
representation is given by 

p{t) ^ dfi{af)dfj.{a'j:)p{a},a'f,t)\af){a'j:\, (23) 



where dfi{a) = ^%|^e I"' is the integrate measure in 

the complex space of the coherent state \a) = e"° |0) [28| . 
and p(af , a'r, t) = {af\p{t)\a'r) can be obtained from the 
exact master equation: 



p{a*f,a'j,t) =A{t) / dfi{ai)dp{a'f)p{a*,a'^,to) 

X exp {a*jB{t)ai + a*jC(t)a'j: 

+ a*D{t)a, + a'*B*{t)a'f}, (24) 

where p(a*,a^,io) = {o-i\p{to)\<^i) is the initial state p(io) 
in the coherent state representation. The exponential 
function in the integral is the propagating function of 
the reduced density operator, which fully determines the 
time evolution of the reduced density matrix in the coher- 
ent state representation. The time-dependent coefficients 
A{t), B{t), C{t) and D{t) in the propagating function are 
listed in table |T] for various cases. Thus, for a given ini- 
tial state p{to), the corresponding analytical solution of 
the reduced density operator for an arbitrary coupling 
to the reservoir at an arbitrary temperature can be ob- 
tained from Eqs. (|23p and (|24l) . Here we shall consider a 
few typical initial states: the vacuum state, the coherent 
state and the mixed state. 



1. Initial vacuum state 

First, we consider the cavity initially in the vacuum 

state, 



p(io) = |0)(0|. 



(25) 



From Eqs. (gSl) and ([H]), it is not difficult to obtain the 
exact reduced density operator at arbitrary time t, 



p{t)^A{t)Y^[C{t)r\n){r 



n=t) 



E 



WT 



-[l + i;(i)]"+ 



Y\n){n\ 



(26) 



This is a completely mixed state in terms of the Fock 
space of the cavity, with the average photon number 



n{t) ^v{t). 



(27) 



That is, there is no existence of a coherent field in the 
cavity if it is initially empty and also no external driv- 
ing field is applied to the cavity. It indicates that the 
coupling to the reservoir makes the cavity continuously 
but randomly gain the energy (photons) from the reser- 
voir, and eventually reach to a steady mixed state with a 
constant average photon number n(t) = v{t — >■ cx)). How- 
ever, we must point out that this steady limit is generally 
different from the BM limit since v{t -^ cxd) y^ n{u}o,T). 
The difference is a manifestation of the non-Markovian 
effect that we shall demonstrate numerically in the next 
section. 

At zero-temperature limit, v(t) — 0, then 



Pit) 



T=0 



p{h)- 



(28) 



In other words, the vacuum state remains unchanged if 
the reservoir temperature is zero. This is a trivial physi- 
cal consequence since at zero temperature, the reservoir 
is also in the vacuum state. Then no photon can be ex- 
changed between the cavity and the reservoir so that the 
cavity remains unchanged. On the other hand, in the 
BM limit, we have 



PBM{t) = X! 



[vBM{tT 



[1 

E 



fBM(i)]"+ 

[n(a;o,r)] 



[l + n{ujQ,T)] 



— rTl'^)(^l 

j\n){n\, 



n+l 



(29) 



which is a thermal equilibrium state with the average 
photon number nsMit) — >■ rL{ujQ,T) at temperature T. 
It is expected that the solution of Eq. (|26|) can be very 
different from that of Eq. (1291) in the strong coupling 
regime. 



2. Initial coherent state 



Next, we consider an initial coherent state, 



Pito) 



^e-|ao|- 



|ao)(ao| 



(30) 



From Eqs. (|23|) and (|24|) . we find that the reduced density 
operator at time t becomes 



^<')--{^} 



..(() 



— n I- 



Ht)] 



n=0 

a{t) 



a{t) 



n,- 



(31) 



where 



"(*) 



,n) =exp[ 



l + v{ty VX'"' l + v{t) 

°^ ' a^]\n) is defined as a gener- 
u{t)a{). It is interest to 



l+i;(t)'"'/ — ^^f^i+v(ty 
aUzed coherent state, and a{t) 
see that Eq. (|3ip is indeed a mixed state of generalized 
coherent states | iV'Jus , n) . The average photon number 
in this state is given by 



i{t) = \a{t)\^ +v{t), 



(32) 



as we expected. 

In the weak coupling regime, u{t) generally decays to 
zero. The steady state limit of the above state will be the 
same as that in the case 1, and the asymptotic average 
photon number is also given by n{t) = v{t — ?> oo). The 
corresponding reduced density operator asymptotically 
becomes a completely mixed state of Eq. (I25|) . This so- 
lution shows an exact decoherence process in the cavity. 
The decoherence arises from two sources, the mean mode 
amplitude damping characterized by the decay behavior 
in terms of the propagating function through the solution 
a{t) = u{t)ao, and the thermal- fluctuation- induced noise 
effect characterized by v{t), as we can see from Eq. (J3T|) . 
The later describes a process of randomly loss or gain 
thermal energy from the reservoir, up to the initial tem- 
perature of the reservoir. 

However, in the strong coupling regime, u{t) may not 
decay to zero, as we shall show explicitly in the numerical 
calculation in the next section. Then the reduced density 
operator remains as a mixed coherent state. On the other 
hand, at zero-temperature limit {v{t) = 0), the reduced 
density operator at time t is given by 



p(<)|T=o-e-l"WI'|a(^))(a(i)l- 



(33) 



In other words, the cavity can remain in a coherent state 
in the zero temperature limit. These two features [u{t) 
may not decay to zero in the strong coupling regime and 
v{t) = at T = 0] indicate that enhancing the cou- 
pling between the cavity and the reservoir and meantime 
lowing the initial temperature of the reservoir can signif- 
icantly reduce the decoherence effect in cavity system. 

On the other hand, in BM limit the reduced density 
operator is always reduced to a thermal state after a long 
time: 



Pbm(^ 



CXD 



[n{^o,TT 



OO 



^\^)inl (34) 



This is because UBmit — >■ oo) must approach to zero, 
see Eq. p^ . It shows that the results in BM limit can 
be qualitatively different from the exact solution. These 
analytical results reveal the underlying mechanism of 
reservoir-induced decoherence in cavity dynamics, which 
may provide some insight how to control the decoherence 
dynamics in such systems. 



3. Initial mixed state 

The last special case we shall consider is the cavity in 
an initially mixed state. 
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with the average initial photon number n(to). Using the 
solution of the master equation, we have 
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where 
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is the average photon number of the cavity during time 
evolution. It shows that if the cavity is initially in a 
mixed state, it will remain in a mixed state, with differ- 
ent particle number distribution varying in time. In the 
weak coupling regime, u{t) approaches to zero so that 
the reduced density matrix at steady limit will reach to 
the same state as in the other two cases. In the strong 
coupling regime, u{t) may not decay to zero so that the 
average photon number in cavity n(i) can be very differ- 
ent from the weak coupling regime, namely very different 
from the BM limit, even though both the exact reduced 
density matrix and its BM limit are mixed states. 

In the next section, we will numerically demonstrate 
these dynamics for some specifically given spectral den- 
sity Jiuj). 



IV. NUMERICAL ANALYSIS OF THE EXACT 
NON-MARKOVIAN DYNAMICS 

To explicitly see the non-Markovian memory effect in 
the cavity system when it strongly couples with a reser- 
voir, we consider a general spectral density of the bosonic 
environment 



J (w) = 2TTr](jj — 
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(38) 
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FIG. 1: Comparison of the exact solution of \u{t)\ for sub-Ohmic (solid Magenta line), Ohmic (dashed red line) and super- 
Ohmic (dash-dotted blue line) with the corresponding BM limit (dotted black line) in the weak coupling {rj < 0.3) and the 
strong coupling {-q > 0.4) regimes. Here we have taken the parameters cjq ~ 21.5GHZ, ujc — (No- 



where r] is the dimensionless coupling strength between 
the system and reservoir, ujc is the cutoff frequency of the 
spectrum. The parameter s classifies the environment as 
sub-Ohmic (0 < s < 1), Ohmic (s = 1), and super-Ohmic 
(s > 1). In the following, we set the value of s to 1/2, 1, 
and 3 respectively. 

For a single-mode cavity with the frequency in the op- 
tical regime, it is well known that the BM approximation 
works very well for a usual thermal environment. How- 
ever, for a structured reservoir, such as the CROW, the 
coupling strength between the cavity and reservoir can be 
controlled by changing the geometrical parameters of the 
defect cavity and the distance between the cavity and the 
CROW 15} . Then the BM approximation or the pertur- 
bative approximation must be reexamined. Furthermore, 
when the cavity frequency lies in the microwave regime, 
the temperature of the environment also becomes non- 
negligible. Specifically we take the cavity frequency to 
be ojQ = 21.5 GHz, i.e., hojo = 13.83^eV, and the tem- 
perature is taken at T = 2K so that ksT = 172.3/ieV 
« 12.5hLLiQ [2g. With these experimental parameters as 
input, we numerically calculate the exact dissipative dy- 
namics of the cavity for three different spectral densities 
with different coupling strength rj in Eq. ([38|) . In addi- 
tion, the cutoff frequency lJc in Eq. psp is taken roughly 
the same order as the cavity frequency wq, i.e. Wc ~ ljq. 
The detailed numerical results are plotted in Figs. [TJIHl 
Note that the BM limits are the same for three differ- 
ent spectral densities when we take the cut-off frequency 
LOc = loq since UBuit), VBuit) and UBuit) only depends 
on K = J{luq)/2 = i:rjWQe~^ for all three different spectral 
densities [see Eq. ^]. From Eqs. (fT9 l) . ([2T1) and (f22|) . 
we can analytically know that u(t), v(t) and n(t) in the 



BM limit are monotonous change in time as the increases 
of the coupling rj. However, in the exact cases, we will 
show that the results have qualitative changes from the 
weak to strong coupling regimes. 

In Fig. [TJ we show the exact solutions for the abso- 
lute value of u{t) (i.e. the amplitude of the propagating 
function u{t) which characterizes the time evolution of 
the amplitude of the cavity field through the relation 
(a(t)) — u{t){a{to))) for Ohmic, sub-Ohmic and super- 
Ohmic spectral densities in the weak and strong coupling 
cases with a comparison to the BM limit. For a given 
coupling strength 77, comparing the behavior of the ex- 
act u{t) of different spectral densities with its BM limit, 
we see that in the very weak coupling limit (r/ = 0.02), 
the difference between the exact amplitude of u(t) and 
the BM limit is very small for all the three spectral den- 
sities. With the increasing of 77, the difference becomes 
more and more visible. The large difference between the 
exact result and the BM limit for a strong coupling is a 
significant manifestation of the non-Markovian memory 
effect. 

For a given spectral density, comparing the behavior 
of \u{t)\ among different coupling strengths, we find that 
roughly for 77 < 0.3, \u(t)\ decays almost monotonously 
for all the three spectral densities, except for a short-time 
oscillation in the beginning. In general, a weak coupling 
to the reservoir always induces an amplitude damping 
to the cavity field. However, when 77 > 0.3, besides a 
short-time oscillation and decay, \u{t)\ may approach to 
a nonzero stationary value. In other words, in the strong 
coupling regime, the non-Markovian memory effect can 
result in a long-time qualitatively change to the time evo- 
lution of the cavity field, namely it changes the cavity 




FIG. 2: Comparison of the exact solution of K(t) for sub-Ohmic (solid Magenta line), Ohmic (dashed red line) and super-Ohmic 
(dash-dotted blue line) with the corresponding BM limit (dotted black line) in the weak coupling {i] < 0.3) and the strong 
coupling {rj > 0.4) regimes. Here we have taken the parameters cjq = 21.5GHZ, Wc = t^o- 



field from a pure damping process in the weak coupling 
regime to a dissipationless process in the strong coupling 
regime. In particular, this qualitative change in the ex- 
act solution of u{t) becomes the most significant for the 
super-Ohmic case, and then for the sub-Ohmic and the 
Ohmic cases, in comparison with the BM limit, as shown 
in Fig. m 

In order to see cleanlier the short- and long-time be- 
haviors of u{t), we plot K{t) = —Re[u{t)/u{t)] in Fig. [2] 
which is the decay coefficient in the master equation ([2]) 
that describes the energy dissipation of the cavity. The 
short-time rapidly increasing of K{t) in the beginning in- 
dicates a fast decay of the cavity field, in agreement with 
the results shown in Fig. [TJ However, with increasing 
the coupling ry, in particular for rj > 0.3, we see that 
the decay coefficient K(t) shows a very different time- 
dependence. It has a large oscillation between an equal 
positive and negative bounded value and then approaches 
to zero. The oscillation indicate that the cavity dissipates 
energy into the reservoir first and gains the energy back 
from the reservoir. Then the zero steady value indicates 
that after the oscillation, the cavity becomes dissipation- 
less. This is why \u(t)\ has a nonzero steady value, as 
we seen in Fig. [1] \u(t)\ of the super-Ohmic reservoir 
maintains a largest nonzero steady value in the strong 
coupling regime because its short-time oscillation occurs 
in the shortest time within which only the smallest en- 
ergy is dissipated. While, in the BM limit, k keeps in a 
nonzero constant, namely the cavity is always in a dissi- 
pation state. 

Therefore, we can conclude that the super-Ohmic 
reservoir contains the strongest non-Markovian memory 



effect shown in the long-time dissipationless process for 
the cavity field, while the sub-Ohmic reservoir involves 
the strongest short-time oscillating non-Markovian dy- 
namics. In the previous investigations on non-Markovian 
effect, one has observed the short-time oscillation be- 
havior in the dissipation processes in the weak coupling 
regime Q. The long-time dissipationless stationary be- 
havior in the strong coupling regime may only be revealed 
from a nonperturbation theory, as we shown here. 

To see clearly the significance of the temperature ef- 
fect of the non-Markovian memory dynamics, we show in 
Fig. [3] the time evolution of the correlation function v{t) 
in the weak and strong coupling regimes and compare 
the results with the BM solution again. The correlation 
function v(t) characterizes the temperature-dependent 
noise effect (the average photon correlation through the 
reservoir). Meanwhile, v{t) is also the main contribu- 
tion to the average photon number inside the cavity in- 
duced mainly by thermal noise dynamics, see Eq. (fT8|) . 
It is shown in Table |I] that, in the zero-temperature limit 
T — >■ 0, v{t) — 0. However, at a finite temperature, it 
will behave very differently for different couplings. 

In the very weak coupling limit (ry = 0.02), the exact 
v{t) is almost the same as the BM solution, except for the 
sub-Ohmic spectral density in which v{t) shows a long- 
time oscillation, as a weak non-Markovian memory effect. 
With the increasing of 77, the quantitative difference be- 
tween the exact solution and its BM limit is enlarged for 
all the three spectral densities. Comparing the behavior 
of v{t) among different 77 for a given spectral density, we 
find that, similar to the behavior of |u(i)|, a short-time 
oscillation (almost invisible in Fig. [3]) of v{t) exists for 
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FIG. 3: Exact solution of v(i) for sub-Ohmic (solid Magenta line), Ohniic (dashed red line) and super-Ohmic (dash-dotted blue 
line) with the corresponding BM limit (dotted black line) in the weak coupling (r\ < 0.3) and the strong coupling (rj > 0.4) 
regimes. Here loq — 21.5GHZ, ojc = too, and T = 2K. 
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FIG. 4: Long time behavior of u{t), v(t) and n(i) for sub-Ohmic spectral density with rj = 0.1. Figure (a) shows Re[u(f)] 
(red solid line) and Im[it(i)] (dashed blue line). Figure (b) shows v(t) (red solid line) and n(t) (dashed blue line). Here 
ojo = 21.5GHZ, ojc = a;o, T = 2K and n(fo) = 50. 



all the three spectral densities. In the long-time limit, 
for different coupling strengths, the stationary values of 
the exact v{t) for all three different spectral densities can 
be very different from the BM result. The numerical re- 
sult in Fig. [3] clearly shows that t;(i — > oo) is far away 
from the BM limit tJEMl^ — > oo) = n{u}o, T) in the strong 
coupling regime, as we discussed analytically in the last 
section. 

In particular, continuously increasing the coupling 
strength r] shows that the sub-Ohmic case manifests a sig- 
nificantly qualitative difference from the BM limit. The 
time oscillation behavior of v{t) for the sub-Ohmic case is 
very strong and maintains for a longer time when the cou- 
pling rj becomes larger. This long-time oscillation behav- 
ior comes actually from the long-time periodical oscilla- 
tion of u{t) , as we shown in Fig. 21 However, in the strong 
coupling regime, v{t) of the super-Ohmic spectral density 



increases slower than the BM limit with increasing rj but 
the overall time evolution is qualitatively the same as in 
the BM limit. In other words, the temperature-induced 
noise effect becomes very important for sub-Ohmic case 
and then for the Ohmic case while it may be minor for 
the super-Ohmic reservoir in the strong coupling regime, 
which is quite different in comparison with the solution 
in the weak coupling cases, see Fig. [31 

To show the totally non-Markovian memory effect dis- 
tributing in the amplitude damping and the noise dy- 
namics, we examine the average photon number inside 
the cavity which is analytically given by Eq. (J18p . From 
equation ([T^. we see that the average photon number 
contains two terms, the first term is the decay of the 
initial average photon number determined by the propa- 
gating function u{t), and the second term is just the cor- 
relation function v{t) as a noise effect which sensitively 
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FIG. 5: Time evolution of the average particle number n{t). Here luq — 21.5GHZ, ujc — ojo, T = 2K, and n(to) ~ 50. 




FIG. 6: Time evolution of the average particle number n{t). Here uo = 21.5GHZ, 



T = 0.002K, and n{to) = 50. 



depends on the initial temperature of the reservoir. The 
numerical solutions based on the exact solution (|18|) for 
three different spectral densities are plotted in Fig.[5]with 
a comparison to the BM result. 

In the weak coupling regime (roughly ry < 0.3), as we 
see the exact solution of n{t) shows a monotonic decay, 
similar to that in the BM limit. However, in the strong 
coupling regime (roughly rj > 0.3), a revival process oc- 
curs in the exact n(t), namely the average photon number 
decays faster in the very beginning and then after a short- 
time oscillation (or no oscillation) , it revives till it reaches 
to a steady value. This behavior does not show up in the 
BM approximation. In fact, the plots in Fig.[5]show that 
the time evolution of the average photon number under- 
goes a critical transition from the weak to strong cou- 



pling regimes. The transition regime here corresponds 
roughly to i] ~ 0.3. This critical transition results in 
a competition between the non-Markovian memory ef- 
fect induced dissipationless phenomena with the thermal 
noise dynamics in the strong coupling in open systems. 

Furthermore, this critical transition can still be seen 
when we reduce the initial temperature of the reservoir 
to a very low value, see Fig. [6] where T — 0.002iir. The 
critical transition is clearly shown up for the super-Ohmic 
case. For the Ohmic and sub-Ohmic cases, the critical 
phenomena is not so significant but the average photon 
number can reach to a very small but nonzero steady 
value in the strong coupling regime. This result indi- 
cates that the initial temperature of the reservoir may 
serve as a sensitive control parameter to control the co- 



11 



herence photon number in cavity dynamics in the very 
low temperature regime, as we have also discussed in the 
analytical solution presented in the last section. 

Put all these analysis together, we find that the non- 
Markovian memory effect can qualitatively change the 
dissipation dynamics of the cavity field in the strong cou- 
pling regime, in particular for the super-Ohmic reservoir. 
Meantime, the non-Markovian memory effects play an 
significant role to the thermal noise dynamics, in par- 
ticular for the sub-Ohmic reservoir. These interesting 
phenomena worth further investigation in other open sys- 
tems. 



V. CONCLUSION 

In summary, we have solved analytically the exact mas- 
ter equation, and obtained the general expression for the 
reduced density operator of the cavity as well as the 
general solution of the mean model amplitude and the 
average photon numbers in the cavity. We take three 
different initial states: the vacuum state, the coherent 
state and the mixed state, to show the different non- 
Markovian time evolution of the cavity state. We find 
that (i) The solution of the exact master equation is very 
different from the BM approximation due to the non- 
Markovian memory effect. In the exact non-Markovian 
case, different initial states may result in different evo- 
lution states and different steady states. While the BM 
solution at steady state limit are the same, independent 
of the initial states, (ii) For the exact non-Markovian evo- 



lution process, temperature effect can play an important 
role, and it drastically changes the thermal noise dynam- 
ics, (iii) The decoherence of cavity coherent field arises 
from both the reservoir-induced damping effect (energy 
dissipation) and the temperature-dependent noise effect. 
Both of them can be controlled by varying the cou- 
pling between the cavity and the reservoir and lowing 
the temperature of the reservoir, as shown in our exact 
solution. These exact analytical cavity dynamics show 
the extreme importance of the non-Markovian effect and 
are numerically demonstrated. Moreover, the numeri- 
cal results show that, the non-Markovian memory effect 
qualitatively changes the amplitude damping behavior 
of the cavity field as well as the thermal noise dynamics 
from the weak coupling regime to strong coupling regime, 
which does not occur in the BM limit. In particular, we 
show that the non-Markovian memory effect leads to a 
long-time dissipationless process to the cavity field in the 
strong coupling regime, and meantime it results in a crit- 
ical transition dynamics when we vary the cavity system 
from the weak to strong couplings with the reservoir. 
Further investigations of all these phenomena in other 
open systems are in progress. 
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